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A DETERMINING FORM FOR THE 2D NAVIER-STOKES 
EQUATIONS - THE FOURIER MODES CASE 



O I CIPRIAN FOIAS 1 , MICHAEL S. JOLLY 2 , ROSTYSLAV KRAVCHENKO 3 '*, 

CN ■ AND EDRISS S. TITI 4 

^ ' Abstract. The determining modes for the two-dimensional incompressible 

Navier-Stokes equations (NSE) are shown to satisfy an ordinary differential 
equation of the form dv/dt = F(v), in the Banach space, X, of all bounded con- 
tinuous functions of the variable s £ M with values in certain finitc-dimcnsional 
linear space. This new evolution ODE, named determining form, induces an 
infinite-dimensional dynamical system in the space X which is noteworthy for 
■ two reasons. One is that F is globally Lipschitz from X into itself. The other 

' is that the long-term dynamics of the determining form contains that of the 

' NSE; the traveling wave solutions of the determining form, i.e., those of the 

(-H , form v(t,s) = vo(t + s), correspond exactly to initial data vq that are pro- 

+a ■ jections of solutions of the global attractor of the NSE onto the determining 

modes. The determining form is also shown to be dissipative; an estimate for 
the radius of an absorbing ball is derived in terms of the number of deter- 
mining modes and the Grashof number (a dimensionless physical parameter). 
Finally, a unified approach is outlined for an ODE satisfied by a variety of 
other determining parameters such as nodal values, finite volumes, and finite 
^ ' elements. 
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O i 1- Introduction 



The notion of determining modes was introduced in [5, 16, 17, 24] as a way to 
gauge the number of degrees of freedom for the Navier-Stokes equations (NSE). 
Determining modes can be described as a (finite) subset of Fourier modes, such 
^ ■ that for any solution {u(t) : i e 1} in the global attractor A its projection into 

those modes, denoted by {v(t) : t £ K}, determines the complete solution. In 
particular, the projection of {u(t) : t £ R} into the higher modes, denoted by 
{w(t) : t £ R}, is determined uniquely by {v(t) : t £ R}. Thus, there is a well defined 
map W : v(-) i-> w(-), which will play a major role in this paper. The number 
of determining modes is associated with the spatial complexity of the long-time 
dynamics, and hence suggests the resolution needed for numerical computations 
and data assimilations (see, e.g., [23, 25]). In this sense, estimates for number of 
determining modes is used much like those for the dimension of the global attractor 
[4, 7, 27]. 

If an inertial manifold exists, the high wave-number spatial modes are enslaved 
by a finite set of low wave- number modes [4, 6, 18, 19, 27]. As a consequence, if 
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the low wave-number modes are specified at just one point in time, the high wave- 
number modes are determined for all time. In the case of an inertial manifold, the 
low modes satisfy an ordinary differential equation called the inertial form whose 
phase space is finite-dimensional. The existence of an inertial manifold for the NSE, 
however, remains open. 

Since a finite set of modes, if specified for all time, can determine the solutions 
to the NSE, a natural question is whether they satisfy an ordinary differential 
equation in a phase space of trajectories, the dynamics of which is consistent with 
those of the NSE. In other words, can one find an analog of the inertial form for 
determining modes? We show in Section 5 that the determining modes v for the 
NSE indeed satisfy a determining form, which is an ordinary differential equation 
of the form dv/dt = F(v), where F : X — > X is globally Lipschitz map, and X is 
the Banach space of all bounded continuous functions from R into the linear finite- 
dimensional space of the low wave- number modes, given in (4.6). The key step is 
to extend beyond the global attractor the enslavement of the high modes mapping 
W : v(-) M- w(-), defined above, to the space X given in (4.6). This is done in Section 
4 by analyzing a differential equation satisfied by w(-), and establishing the Lipschitz 
property of W . We then show in Section 6 that the resulting determining form, 
dv/dt = F(v), is dissipative, and estimate the radius of its absorbing ball in terms 
of the number of determining modes, and the Grashof number (a dimensionless 
physical parameter given in (2.26)). 

Though the number of determining modes is finite, the ODE they satisfy, i.e. 
the determining form, induces an infinite-dimensional dynamical system, since it is 
an evolution ODE in the infinite-dimensional Banach space X. Just as the set of 
partial differential equations that comprise the incompressible NSE can be written 
as an evolution equation, involving unbounded operators, in a Hilbert space H 
of functions in the spatial variable, the determining modes evolve as solutions to 
an ODE in a Banach space X of functions in an auxiliary time variable s. A 
significant difference, however, is that the vector field that governs the evolution in 
the determining form ODE is a globally Lipschitz map from the space X into itself, 
whereas in the case of the NSE evolution equation it involves unbounded operators. 
Rather than being a dimension reduction of the NSE, the determining form is in 
fact an embedding of the original dynamics in a larger phase space. In the phase 
space X, the solutions at any moment in time i, represent complete trajectories 
v(t, s), — oo < s < oo. 

For certain initial trajectories the behavior of solutions to the determining form 
is simple to describe. We show in Section 7 that traveling wave solutions to the de- 
termining form can be characterized as those satisfying u(i, ,s) = p(t + s; uq), — oo < 
s < oo, where p(s;uo) denotes the projection into the determining modes of the 
solution at time s to the NSE with initial condition uo £ A. As shown in Section 5 
any solution with an initial trajectory u(0, s), — oo < s < oo which is periodic in s, 
with period T, remains periodic in s, with period T, regardless of whether this initial 
trajectory is on the global attractor of the NSE. As a special case, any solution with 
an initial trajectory independent of s remains so; if v(0,s) = c(0), — oo < s < oo, 
then v(t, s) = c(t), — oo < s < oo. 

We also compare the general long time behavior of the determining form with 
that of the NSE. In Section 9 we consider the sequences v n (t, s) = v(t+t ni s—(t+t n )) 
and w n (t, s) — w(t + t n , s — (t + t n )). We show that v n + w n converges weakly to 
a solution of the NSE with an additional Reynolds force which has only finitely 



A DETERMINING FORM FOR THE NSE 



3 



many modes. The stationary solutions of the determining form satisfy a differential 
algebraic system: an algebraic relation for v, coupled with a differential equation 
for w. Analogs of the energy and enstrophy balances are used in Section 10 to 
show that, as in the case of the NSE, a certain extremal condition on norms of 
the stationary solution occur only for forces that are eigenfunctions of the Stokes 
operator. 

Over the years it was proven that there are several other finite systems of param- 
eters, different than the modes considered here, which also determine the long-time 
behavior of the solutions to the Navier-Stokes equations (e.g., determining nodes, 
finite volumes averages, finite elements, etc. [3, 15, 21, 20, 24]). It is worth stressing 
that the determining form studied here uses only determining modes, on which the 
partial differential operators in the NSE act in an obvious manner; this action is 
seemingly quite involved for other systems of determining parameters. One also 
faces a similar problem in representing inertial manifolds in terms of determining 
parameters other than the Fourier modes. This challenging general parametrization 
issue of the inertial manifolds has been solved in [22] for the determining nodes case 
and later in [3] for the general interpolant case, including the finite elements case 
etc. In a forthcoming paper [12] we will introduce a general, unified construction 
of determining forms based on any finite system of determining parameters. This 
unified approach has also been used in [1, 2] in the context of continuous data 
assimilations and finite-dimensional feedback control, respectively. In section 11 
of this paper we briefly illustrate this general construction for the finite system of 
determining modes considered in this paper. 

2. Preliminaries on the Navier-Stokes equations 
We write the incompressible Navier-Stokes equations 
On 

(2.1) — - vAu+ (u ■ V)m + Vp = $, 

at 

(2.2) divu = 

(2.3) udx = , I $dx = 

Jn Jn 

(2.4) u(t ,x) = u (x) 



subject to periodic boundary conditions with basic domain fl = [0,L] 2 as a differ- 
ential equation in a certain Hilbcrt space H (see [4] or [27]), 

ju{t) + uAu{t) + B{u{t), u(f)) = /, 

u{t) G H , t > , and u(0) = u . 

The phase space H is the closure in L 2 (f2) 2 of the set of all R 2 -valued trigonometric 
polynomials <f> such that 



V • (j> = 0, and / (j){x)dx = 0. 

Jn 

The bilinear operator B and force / arc defined as 

(2.6) B(u,v) =V((u- V)«) , f = V§, 

where V is the Helmholtz-Leray orthogonal projector of L 2 (f2) 2 onto H, and where 
u, v are sufficiently smooth. When the input to B is lengthy, and repeated, for 
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convenience we will use the shorthand notation 
(2.7) B(u) = B(u,u). 

Unless specified otherwise, / G H. The scalar product in H is taken to be 
(u,v) = / u{x) ■ v(x)dx, where a ■ b = Oi&i + 02^21 



with associated L 2 -norm 

1/2 

(2.8) |u| = {u,u) 1/2 = ( / |u(aj)|da; 

where, for a vector a € C 2 , we denote its length by \a\ = (a ■ a*) 1 / 2 , with a* = 
(01,02)- It will be clear from the context when | • | refers to the length of a vector 
in C 2 , as opposed to an L 2 -norm. The operator A = —A is self-adjoint, and its 
eigenvalues are of the form 

27r\ 2 

— J k-k, where fc e Z 2 \ {0}. 

We denote these eigenvalues by 

< A = (2tt/L) 2 < Ai < A 2 < ■ • • 

arranged in increasing order and counted according to their multiplicities, and write 
Wo, W\, W2, ■ ■ -, for the corresponding normalized eigenvectors (i.e., \wA = 1 and 
Awj = XjWj for j = 0, 1, 2, . . .). 

The positive roots of the operator A are defined by linearity from 

A a Wj = Xfwj , for j = 0, 1, 2, . . . 

on the domain 

00 

D{A a ) = {ueH :^A 2q (u, Wj ) 2 < 00}, 
3=0 

for a > 0. 

We take the natural norm on the space V := D(A 1 / 2 ) to be 

1/2 

/ p £ rj rj ^ 

(2.9) ||u|| = \A 1/2 u\ 



^ f 2 d d \ f 00 



Since the boundary conditions are periodic, we can express any element in H as a 
Fourier series 

(2.10) u(x) = ]T u k e iK ° k - x , 

fcez 2 

where 

(2.11) k = A /2 = , wo = , Mfe = u_ fc , 

and due to incompressibility, k ■ Uk — 0. We associate to each term in (2.10) a wave 
number Ko\k\. Parscval's identity reads as 

M 2 = L 2 £o fc .^ fe = L 2 ^K| 2 , 

feez 2 kez 2 
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as well as 

(u,v) = L 2 ^ u k ■ v^ k , 

for v = Vke lK ° k ' x . We define the orthogonal projectors Pn,Qn, Pn ■ H — > 
span{iOj|Aj < (k N) 2 }, by 

Pi \ * - i/ink-X 

\k\<N 

where u has the expansion in (2.10), and Qn = I — Pn- 

Recall the orthogonality relations of the bilinear term (see, e.g., [4, 27]) 

(2.12) (B(u, v),w) = -{B(u, w), v) , u, v, w e V , 
where (-, •) denotes the the action of an element in V' , as well as 

(2.13) (B(u,u), Au) = , ueD(A). 
The strengthened form of the enstrophy invariance 

(2.14) (B(Av,v),u) = (B(u,v),Av) 7 ueV 7 v £ D(A) 
is proved in [9]. Relation (2.14) together with (2.12) implies 

(2.15) (B(v,v),Au) + (B(v,u),Av) + (B(u,v),Av) =0, u,v £ D{A) , 

(2.16) (B(Av,v),u) - (B(v,Av),u) + (B(v,v),Au) = , u,v £ D(A) , 
see [13]. We will also use Agmon's inequality 

(2.17) M» < c A M 1/a |4«| 1/a , u€D(A), 
and one of its consequences 

(2.18) \(B(u,v),w)\ ^caImIIIwII 172 !^ 3 / 2 ?;! 1 / 2 ^! , u £ H, v £ D(A 3 / 2 ), w £ H . 
(See [4, 26, 10]). 

In this paper we will assume / £ D(A>' 2 ) for some specified j £ {0,1,2,3}, 
which is sufficient to guarantee that 

1 (J, l9 ,du , 1 d „ ,,, ,du , . 

in the sense of distributions (see Lemma 1.2 in Chapter 3 in [26] and its proof). 
Taking the scalar product of (2.5) with u, respectively Au, and applying (2.12), 
(2.13), and (2.19) gives the energy, enstrophy balance equations 

(2.20) ~N a + HMI 2 = (/,«), 

(2.21) ~|| U || a + I /|Au| 2 = (/,^«). 
In the scientific literature, 

Kq|u| 2 = 2 times the total energy per unit mass , 

and 

K o!M| 2 = the total enstrophy per unit mass . 
Straightforward applications of the Cauchy-Schwarz, Young, and Gronwall in- 
equalities to (2.21) gives 

(2.22) ||u(t)|| 2 <e-^ t ||u(0)|| 2 + -^(l-e- ,yK «*) . 
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If |/| = 0, we have by (2.22) that ||M(t)|| — >• as t — >• oo. Therefore we assume 
throughout the paper that 

(2.23) |/| > . 

It is clear from (2.22) that the closed ball B of radius 2vkoG, centered at 0, is 
absorbing and invariant, meaning that for any uq there exists T = T(||uo||) such 
that 

(2.24) ||u(t;uo)|| < 2vk G , for all* > T , 
and 

(2.25) \\u(t;u )\\<2vKoG , for aU t > , u G B 
where 

(2.26) G = 44 

is the generalized Grashof number introduced in [16]. It is well known (see e.g. [13]) 
that for the flow to be turbulent, we must have G> 1. 

While most of our results apply for all G > 0, wc will eventually assume G> 1. 
We also have from (2.22) that 

(2.27) 1 1 -SLA 1 1 < vkqG for all u e A , 

where A is the global attractor, i.e. the minimal compact set in V which attracts 
uniformly all solutions of the Navier-Stokes equations starting from the ball {h € 
V : \\h\\ < 2Gvkq}. Recall that this set can also be characterized as 

(2.28) A = {u eH:3 u(t, u ) = solution V t £ E, sup ||u(i)|| < oo}. 

i 

We often denote the solution as simply u(t) when the choice of uq is irrelevant to 
our considerations. It is easily shown using (2.22) that T, the time of absorption in 
(2.24), can be taken to be 

(2.29) nND^^bgWL}. 

We will use the following bounds on the nonlinear term (see [8, 15, 28] and the 
references therein) 

(2.30) \{B(w,u),v)\ < c T \\w\\\\u\\ (in 4 



(2.31) \(B(w,u),v)\ < c B \\w\\\\u\\ In -L—1 \ v \, 

K o\\ w \\ 



K \V 



e\Aw\ 



1/2 



e|,4u| 



1/2 



(2.32) \(B(w,u),v)\ < ct||w||||u|| ( In— ^— j 1 \v\, 

V K o\\u\\J 

(2.33) \(B(u,v),w)\ < c L \unu\\V 2 \\v\\\wnw\\V 2 , 

(2.34) \(B(u,v),w)\ < c L |ur /2 ||u|| 1/2 ||w|| 1/2 |^r /2 |w|, 

(2.35) \(B(u,v),w)\ < c A |«r /2 |^r /2 |k||k| 

(2.36) \{B(u,v),w)\ < c A |w||l^|| |^| 1/2 |^L^| 1/2 
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which are valid whenever the norms involved make sense. One can use these in- 
equalities to extend the domain of B to larger functional spaces. For instance from 
(2.35) one can infer that B(u, v) can be extended by continuity to a map from 
D(A) x V to H . We will denote these extensions also by B. 

3. Specific preliminaries 

We will use the following elementary estimates in several places. 

Lemma 3.1. Let x(t), y(t) be nonnegative functions defined on the interval I, 
where I = [to,t\), for to finite, or (— oo,ti). Suppose that x 2 (t) is differentiable on 
(to,ti), and satisfies 

(3.1) j t x 2 (t) < ~ay 2 (t) + by(t), x(t) < ey(t) 

for t 6 / where a, b, e are positive constants. Then the following hold: 
(i) If x(t) > eb/a for all t € [to,t±), then 

(3.2) x(t) -e-< e -o(*-*o)/2e 2 f x f to ) _ e _); 

a a 

(ii) if x(t) < eb/a, for some t G I, then x(t) < eb/a for all t £ [t,t\); 
(Hi) if x(t) is bounded on (~oo,ti), then x(t) < eb/a for all t £ (— oo,ti). 

Proof. To prove (i) define c(t) := e _1 a;(t) — b/a. Since 

(3.3) y(t)>^- = c{t)+ b -> h - 

e a a 

we have 

dc 2e 2 , , . „ 9 dc n 9 b dc 9 

ac + b ) = 2e 2 c— + 2e 2 -— = -x 2 < - ay 2 + by 

at a at a dt at 

< - a(c + -) 2 + b(c + -) = -c(ac + b). 
a a 

Thus, since ac + b = ax/e > b > 0, 

(3.4) e atl2e \{t) < e ato/2£2 c{t ) 
and (i) follows. 

To prove (ii), suppose that x(t') > eb/a for some t' > t. It follows that both 
-^.x 2 > and x > eb/a hold at some time between t and t' . We have at that time 
by (3.3) 

(3.5) < j t x 2 < - ay 2 + by < 0, 
a contradiction. 

(iii) follows from (i). □ 

Wc also use the following immediate corollary of Lemma 3.1 for the case in which 
the coefficient b varies with time. 

Corollary 3.2. Suppose x(t), y(t) are positive functions defined on the interval I, 
where I = [to,ti) for t finite, or (— oo,t\), and e > is such that x(t) < ey(t) and 

(3-6) j t x 2 {t) < -ay 2 (t) + b(t)y(t) 

for t £ I, with a, b(t) positive, then 
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(i) if x(t) < esup tl< fb(t')/a, for some i £ I, then x(t) < £sup t , <t b(t')/a for 
all t e [t,h) 

(ii) if x(t) is bounded on (—00, t\), then x(t) < eb/a for all t G (— 00, tx). 

Definition 3.3. The projection P^ is called determining if (and only if) for any 
two solutions u\{i),u 2 {t) (t € R) in A we have 

(3.7) P N ui(t) = P N u 2 (t) =^ «i(t) = u 2 {t) 

This is equivalent to the original definition in which the high mode trajectories 
are determined by the limit as t ~ > 00 of their projections [14, 15, 17]. Using the 
original definition, it is shown in [24] that for P/v to be determining, it is sufficient to 
take N — 0(G). For completeness, we now reproduce that estimate using Definition 
3.3. We recall (see [15], Ch.3 Sec.l) that u = u\ — u 2 satisfies 

(3.8) \j t \ u \ 2 + = -(B(u,ui),u) < clIuHIuIHIuiII < ct|«|||u||Gi/«o, 
whence (since Pnu = and so n n N\u\ < \\u\\) 

(3.9) | H 2 + 2HM |2< 2CL ^| M |2 
d 1 \2 o..Ar2\ (, Gc h \ 2 



(3.10) _ N 2 + 2 ^ Ao M__^j|^< 
so that 

(3.11) \u(t)\ 2 < | u ( to )|a e -a^AoCi-G OL /Jv)(t-t ) j v to < t 
Thus if 

(3.12) N > c L G, 

we conclude (by letting to — > —00) that u(t) = 0; that is if (3.12) holds, then Pjy is 
determining. Note that (3.12) implies that if u\(t) = u 2 (t) only for all t < T then 
still u\{t) = u 2 (t). In particular Qnu(-) is uniquely determined by Pnu(-) in case 

(3.12) holds. Then by denoting 

(3.13) q(s) = W(s; P N u(-)) = Q N u{s) (s e M), 

we have that on the set {Pnu(-) : u(-) is a solution of (2.5) in .4} the map 

(3.14) P N u{-) i-> q(s) = W(s; P N u{-)) (s e R) 
is well-defined by the following properties 

(3.15) + uA ^ s ) = Q N f _ Q N B(P N u{s) + q(s)), Vsel 
as 

(3.16) sup ||q(.s)|| < 00, 
where we used the shorthand notation in (2.7). 
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4. Extending the domain of W to C&(R, PnV) 

To ultimately write a differential equation for the low mode trajectories, we need 
to extend the domain of W beyond the projection of the global attractor into PnV . 
We define 

(4.1) p(t,s) = P N u(s + t), (Vs,i6K), 

and apply the projection Pjv to (2.5) to obtain 

(4-2) %-p(t, s) = -vAp(t, s) + P N f - P N B(p(t, s) + W(t + s; P N u{-))). 

at 

But it is clear that 

(4.3) W(t + s; P N u(-)) = W{s; r t P N u{-)) = W(s; p{t, ■)), 
where r t is the shift by t: 

(4.4) T t u(s) = u(s + t) , for all s G ffi . 
Thus (4.2) becomes 

(4.5) ^p(t, s) = -isAp(t, s) + P N f - P N B(p(t, s) + W(s; p(t, •))). 

We now consider the Banach spaces 

X = C b (R, P N H) = {»:Rh> P n H, v(s) continuous Vsel, \\v\\ x < oo}, 

(4.6) 

Y = C b (R,Q N H) = {w : K ^ Q N H,w(s) continuous V s G E, \w\ Y < oo} 

where 

(4.7) \\v\\ x =sup||u(s)|| , \w\ Y =sup|tu(s)| . 

s s 

Let p(t, s) in (4.1) be written as v(t) G X, so that (4.5) can be viewed as the 
"differential equation" 

(4.8) ±v(t) = -vAv(t)+P N f -P N B(v(t) + W(v(t))) (VtGR), 
at 

where throughout v(t) denotes the function s i— > v(t,s) in X and VF(w(t)) denotes 
the function s i-> M^(s;w(t)). We should stress that dv/dt in (4.8) is the derivative 
of a trajectory in the space X. It is also useful to note that, in particular, the 
following partial differential equation holds 

|U(i, s) = P N f - vAv{t, s) - P N B(v(t, s) + W(s; v(t, •))) , 

for every s G K, where d/dt is taken in H . Note that (4.8) is defined only for those 
v{t) that come from p(t, s), i.e. from u G A. 

Now let y(£) be a continuous, piecewise linear function, that is equal to 1 for 
|£| < 2G, for |£| > 3G, where G is the Grashof number given in (2.26), and 
consider 

(4.9) - v (t) + vAv{t) + P N B (<pv(t) + W(v{t))) = P N f , 
where 

(4.10) (p = <p 



Ht)\\ 



We will show that W(v) can be extended to a Lipschitz map X — > Y in such a way 
that (4.9) defines an ordinary differential equation in X, which coincides with (4.8) 
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on the ball of radius 2vk$G in X. By this we mean that v(t) in (4.9) is the solution 
of an equation 

dv 

(4.11) di =F( -^ ( V * eR > veX )' 

where F : X — > X is a locally Lipschitz function. In fact we will show that F is 
globally Lipschitz 

(4.12) ||F(wi) - F(v 2 )\\ x < L F \\vt - v 2 \\x, 

where the constant Lp is independent of v\ , v 2 . This will provide a proof that equa- 
tion (2.5) on the attractor A is transformed into an ordinary differential equation 
in the space X defined in (4.6). 

Let h £ QnH, v be an element of X. We will show that 



(4.13) 



— w + vAw + QjsrBUpv + w, tpv + w) = h 
as 



has a unique bounded solution w(s) £ Y. Consider a Galcrkin approximation for 
77, > N, of equation (4.13), namely 

-w n + vAw n + P n QNB(tpv + W n ,fV + w n ) = P n h 



(4.14) 



ds 



w n (s) = P n w n (s) for s> s , 

with the initial condition w n (so) = 0. Since (4.14) is an ordinary differential equa- 
tion, it has a unique solution w n (s) for s in some maximal interval [so, si), on which 
||iu(s)|| < oo. 

Taking the scalar product of (4.14) with Aw n and using relation (2.15) we obtain: 

(4.15) ~\\w n \\ 2 + v\Aw n \ 2 + tp 2 (B(v,v),Aw n ) - tp(B(w n ,w n ),Av) = (h,Aw n ) . 
2 as 

Note that it follows from the definition of <p that 
HI 

l/KO, 

Using (2.32) and (4.16), we have 

<p 2 \{B(v,v),Aw n )\ < <p 2 c T \\v\\ 2 \Aw n \{\neN) 1 / 2 < 9c t p 2 k 2 G 2 (In eN)V 2 \Aw n \ , 
while using (2.33), ||tu n || < (1/ K N)\Aw n \, \A 3 ^ 2 v\ < (k N) 2 \\v\\ and (4.16), we get 
ip\(B(w n ,w n ),Av)\ = <p\(B(w n ,Av),w n ))\ < c L |w„| \\w n \ \\A 3 ' 2 v\ < 



(4.16) 



HI < 3^k G VtjgX, Va>0. 



Thus 



(4.17) 



G 

3c h v—\Aw n \ 2 . 



1 d 



v\Aw n \ 



2ds n " 

< \h\\Aw n \ + ip 2 \{B(v,v),Aw n )\ +<p\(B(w n ,w n ),Av)\ 



G 

< 3c L ^— I^Tfnl 2 



^v 2 K 2 G + %c T v 2 K 2 G 2 {\neN) 1 / 2 



\Aw n \ 



Now we have by (ii) of Lemma 3.1 (with e = 1/(kqN)), that if N and G are such 
that 



(4.18) 



TV > 3c L G, 
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then 

,,,,, 9o T G 2 (lneW)" 2 + K{G 

<««> IkMII < «. — n _ 3clG m - 

In the following we assume that (4.18) is satisfied, and introduce the notation 

9c T G(lneA0 1/2 + U 

(4.20) K g,n)= n 3clG b . 

Thus for the solution w n (s) of (4.14) with initial condition w n {so) = 

(4.21) IK0)|| <vk GP(G,N) . 

It follows from the general theory of ordinary differential equations that w n can 
be extended to the solution of (4.14) defined on [so,oo). Therefore the classical 
proofs of existence for solutions of equation (2.5) (see e.g. [4] Ch.9 and Lemmas 8.2 
and 8.4 in Ch.8) apply almost verbatim to 4.13. It follows that (4.13) has a unique 
solution on [so,oo) satisfying w(sq) = 0. Now denoting by w^ n '(-) the solution of 
(4.13) on [— n,oo) satisfying w^(— n) = 0, we have 

(4.22) \\w {n) (s)\\ < K Q vGp{G,N) for all s £ [-n,oo). 

It is also easy to show using (4.17) that for any fixed, bounded interval [so,Si] we 
have 

(4.23) sup / 1 \ Aw (n) {s)\ 2 ds < oo. 

n J s 

It is well known how to prove that there exists a subsequence {rij}°^L 1 such that 
w ( n j) converges to a solution w(s) of (4.13) in L 2 ([so, si]; V) satisfying 

(4.24) ||tu(s)|| < K isGf3(G,N) for all s e [s ,Si] , 
and moreover, there exists C = C(G), such that for si — sq < 1, 



(4.25) 



/•S2 

/ |Aw(s)| 2 < Cvk^vG for all si,S2 such that s\ — sq < 1 

J si 



(see again [4]). 

We note that if w(s) is any bounded solution of (4.13), defined on (— oo, si] then 
in exactly the same way we get the analog of (4.17) for w, from which by (iii) of 
Lemma 3.1 we get the bound 

(4.26) \\w(s)\\<vk G/3(G,N) for all s G (-oo, si]. 

We now prove that this solution, which is bounded in the || • || norm, is unique 
for sufficiently large N. Let wi,W2 be two bounded solutions of (4.13). Denote 
w = wi — u>2 ■ Then 

— w + vAw + Q N B(ipv + wi,(pv + wi) - Q N B(ipv + w 2 ,<pv + w 2 ) = 0. 
as 

Taking the scalar product of the above equation with w we get 

~M 2 + v\W\ 2 

= —(B((pv + wi,(pv + wi),w) + (B(ipv + w 2l ipv + W2),w) 

= —{B(w, ipv + u>i), w) — (B(ipv + u>2, w), w) = —(B(w, ipv + wi), w). 
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We apply (2.33), (4.16), and (4.24) to obtain 

\(B(w,<pv + Wl ),w)\ < clMIMIMMI + IKH) < G(3 + p&NVdKovWWwW 
so that 

(4.27) ~\w\ 2 + v\\w\\ 2 <G(3 + j3(G,N))c h K v\\w\\\wl 

2 as 

Now we use 

(4-28) \w\ < -^\\w\\ 

k N 

to get 



d l..,2 , o„2.. Ar2 n „ „3 + p(G,N) s ,, 2 



(4.29) — \w\ 2 + 2k 2 q vN 2 (1 - c h G ^ ' ' )\w\ 2 < 0. 

Since \w\ is bounded on (—00,00), we find that w = 0, when 

(4.30) N > clG(3 + /3(G, N)) . 

We next show that w(s) is continuous in s. Integrating (4.13), we have 

\w(sz) — w(si)\ — I / h — QnB(v>v(s) + w(s)) — vAw{s) ds\ . 

J Si 

Using (2.31), (2.32), (2.35), together with (4.16) and (4.24), we find that 

\B((pv)\ < 9caKo) 2 G 2 , \B(tpv, w)\ < 3c B (^o) 2 G 2 /?(G, N) , 
\B(w,<pv)\ < 3c T (v K0 ) 2 G 2 p(G,N) , \B(w)\ < Sffi(yG/3(G, N)) 3 / 2 \Aw\^ 2 
It follows that (for < S2 — si < 1) we have 

\w{s 2 ) - w{ Sl )\ <{\h\+ [(9ca + 3I3(G,N)(c b +ct)G 2 +C(G)] N 2 }( S2 -si) 
+ ] ^4 /2 - 7/4 (G/3(G,iV)) 3 / 2 [G(G)] 1 /4 (s2 _ Sl) 3/4 . 

We work with both the | • |- and j| • ||-norms to allow for h (hence /) to be in H, 
rather than in V. 

Note that (4.30) implies (4.18). This concludes the proof of the following lemma. 

Lemma 4.1. Assume that N satisfies (4.30), where f3(G,N) is defined in (4.20). 
Then the equation (4.13) has a unique weak solution defined on the whole real line, 
which is bounded in the || • || norm. Moreover this weak solution is in Y with 

(4.31) \\w(s)\\<uk G(3(G,N) VseR, and \w\ Y < vG ^°^ N \ 

This lemma allows us to define the (nonlinear) operator W which takes v & X to 
the unique solution of the corresponding equation (4.13). By Lemma 4.1 we have 
that 

(4.32) W{s;T t (v{-))) =r t W(s; <)). 

It is also useful to have a bound for \Aw\ if \\h\\ is finite. Consider the following 
a priori estimates of Aw. Taking the scalar product of (4.13) with A 2 w we have 

(4.33) -^-\Aw\ 2 + v\A 3 / 2 w\ 2 + {{B(<pv + w),A 2 w) = [h 1 A 2 w) < \\h\\\A 3 / 2 w\. 
2 ds 
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Using (2.16), expanding, then applying (2.12) (three times) and (2.14), we obtain 

{{B{<pv + w, ipv + w),A 2 w) 

= (p 2 {B(v,Av),Aw) + (p(B(v,Aw),Aw) + tp(B(w,Av),Aw) + (B(w,Aw),Aw) 

- ip 2 (B(Av, v), Aw) - ip(B(Av, w), Aw) - ip(B(Aw, v),Aw) - (B(Aw, w), Aw) 
= -tp 2 {B(v, Aw),Av) - ip{B(w, Aw), Av) 

- ip 2 (B(Aw, v),Av) - (p(B(Av, w), Aw) - (p(B(Aw, v),Aw) - (B(Aw, w), Aw) 

so that by applying for the first term (2.30) and for the rest of terms (2.33), followed 
by (4.16) and (4.31), we have 

\((B(ipv + w,ipv + w), A 2 w)\ 

< CT^IIulll^wllAuKlne^^ + CL^Iwl^ll^ll^i^^ii^ii/si^^ii/a 

+ c hV 2 \Aw\ 1 / 2 \A 3 ' 2 w\ l ' 2 \\v\\\Av\ 1 / 2 \A 3 l 2 v\ 1 ' 2 
+ c L ^\Av\ 1 / 2 \A 3 / 2 v\^ 2 \\w\\\Aw\ 1 / 2 \A 3 / 2 w\ 1 / 2 
+ c L v3|Aw||^ 3/2 u;|||w|| + c L |A 3/2 u>||,4u;|||u>|| 

< ctMMI) 2 K N(lneN) 1/2 \A 3/2 w\ + c h <p\\v\\ 



\A 3 I 

+ clMM|) 2 koN|^ 3/2 H + CL<p\\v\y- — 



k N 



2 w\ 2 



k N 

\A 3 ' 2 w\ 2 Q .„ ^JA 3 / 2 w\ 2 



+ c L tp\\v\y- -i- +c l k is(3(G,N)G 



KqN ' ' K N 

<c L [9 + [3(G, N)]v^\A 3 / 2 w\ 2 + 9^ 2 K 3 G 2 iV(c T (lne^) 1 / 2 + c L )\A 3 / 2 w\ 
Thus if 

(4.34) N > c L G(9 + (3(G,N)) , 

we can apply Lemma 3.1 (iii) to obtain 

\\h\\ + 9is 2 K 3 G 2 N(c T (liieN) 1 / 2 + CL ) 



(4.35) \Aw\ < 



vk {N- CIj G(9 + P(G,N))} 



Using the Galerkin approximation we get the bound (4.35) for w n , provided 
(4.34) holds. Therefore, w n € L°°(0, T; D{A)). From (4.33) we obtain that each 
w n belongs to a bounded set of L 2 (0, T; D(A 3 / 2 )) for any T. Therefore as usual 
choosing a subsequence that converges weakly in L 2 and weakly* in L°° , and taking 
the limit as n — > oo, we obtain the solution of (4.13). Note that (4.34) implies 
(4.30). Therefore the following lemma is true, 

Lemma 4.2. Suppose h e D{A 1 / 2 ). If G and N satisfy (4.34), then the solution 
to (4.13) satisfies 

\\h\\+9v 2 K 3 G 2 N(c T (lneN)V 2 + c L ) 

(4 - 36) lAwl - ^{N-cLGid+mm)} ■ 

We note that the analog of the results in this section with w <G Cb((— oo, T], Q^H) 
instead of Cj,(R, QjyH) in this section remains valid for v € Cj,((— oo, T], P^H). 
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5. The determining form 
5.1. Lipschitz property of W . Let v±, V2 be in X, and w% = W(vi), W2 = W(v2). 
Denote 7 = w\ — w 2 , 8 — v\ — v-i-, i'i = tp {j^rf^j u» 5 i = 1,2, 8 = V\ — v 2 - Subtracting 
equation (4.13) for W2 from that for w\, we obtain 

(5.1) — j + vAj + Q n B(5 + j,v 1 +w 1 ) + Q N B(v 2 + W2,5 + 7) = 0. 
as 

Taking the scalar product with A'y, we have 

(5 2) \i M2 + " m ~ 6 ' il + Wl ^ A ^ + ^ + w 2^),^7)l 

+ |(B(7,t5i),A7)| + \B(v 2 ,l),Aj)\ + \(B(i,wi), A 7 )| + \B(w 2 ,j), A 7 )|. 

To estimate the right hand side we apply (2.31) to the first term, (2.32) to the 
second, (2.34) to the third term, (2.12), (2.14), and (2.36) to the fourth term, and 
(2.35), (2.34) respectively to the fifth and sixth term. Using then (4.16), (4.24) and 
(4.28), we have 

' " 7 || 2 + H^7| 2 < G(3 + P(G,N))k q v(c t + ^fSlWA^lneN) 1 / 2 



(5.3) 2 ds 



|2 G , , .Gp(G,N). A |2 



+ 3v(c A + c L )|^ 7 | 2 - + v(c A + cl) y ' \A 7 \ Z , 



and thus by Lemma 3.2 (ii) we have that 

(5.4) < ""/^Z':? 1 ™ p 

N - G(3 + fi(G,N))(c A + c h ) s >< s 

for G, N such that 

(5.5) N>G(3 + p(G,N))(c A + c L ) . 
From now on we will consider G and N such that 

(5.6) N>G(9 + (3(G,N))(c A +c L ) 

which implies (4.18), (4.30), (5.5) and (4.34). 

We derive from (5.4) the following two corollaries 

Corollary 5.1. Suppose that G and N satisfy (5.6). Then ifvi,v 2 € X, Vi(s) = 
112(3) for all s < so we have that W{s, v\) = W{s, v 2 ) for all s < sq. 

Proof. Indeed, if s < sq, it follows that 8(a) = for all a < s, thus (5.4) gives that 
7 (s)=0. □ 

Corollary 5.2. Suppose that G and N satisfy (5.6). Then W : X — s- Y is a 
Lipschitz function with Lipschitz constant Lw satisfying 

. G(3 + /3(G, iV))(lneiV) 1 / 2 (c r + c B ) 



(5.7) L w < 4- 



iV-G(3 + /3(G,iV))( CA +c L ) 



Proof. If we have both ||t>i|| > 3i/kqG and ||«2|| > 3G, then 8 = 0. If, say ||wi|| < 3G, 
then 

II* - fell < \<p ( M ) - 9 ( M ) HKII + * (M) , K _ V2 || 

< — ^HK'ill - ||w 2 |||3i/koG + ||f ! - u 2 || 
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thus 

(5.8) \\6\\ < i\\6\\ . 

□ 

Theorem 5.3. Suppose G and N satisfy (5.6). The equation (4.9) is an ordinary 
differential equation in the Banach space Cf,(K., PnH). 

Proof. It suffices to show that F(v) = P N f - vAv - P N B(v + W{v), v + W(v)) is 
Lipschitz. To estimate Lp, the Lipschitz constant of F, we apply (2.30), (4.16), 
(4.24) and Corollary 5.2 to obtain for every (nonspecificd) s g R 

n(s) = \P N B(v 1 + W{v 1 ),v 1 + W(v 1 ))-P N B(v 2 + W(v 2 ),v 2 + W{v 2 ))\ 

< \PnB(v! + W{v 1 ),v 1 + W{v t )) -P N B(v 1 + W(v!),v 2 + W{v 2 ))\ 

+ \P N B(v! + Wivt), i 2 + W{v 2 )) - P N B{i 2 + W(v 2 ),v 2 + W(v 2 ))\ 
= \P N B(h + W(«i), vi - v 2 + W{v x ) - W(v 2 ))\ 
+ \P N B(v! ~ v 2 + W{vi) - W(v 2 ),v 2 + W(v 2 ))\ 

< c T (\\ii + W(v!)\\ + \\v 2 + VK(v 2 )||)(lneA0 1/2 ||£i - «a + Wfa) - W(v 2 
It follows that 

sup?/( S ) < 2c t vk G(3 + 0(G, ^))(lne7V) 1 / 2 (3 + L w )\\vi - v 2 \\ x 

and therefore 

WFM-F^Wx < Ko N skp\F( Vi { S )) ~ F{v 2 {s))\ 

< [vk 2 N 2 + 2c T iyn 2 Q NG(3 + /3(G,N))(\neN)^ 2 (3 + Lw)}^ - v 2 \\ x 

□ 

We denote the solution operator for (4.9) by Sdf, i- e - 5 , DF(i)«o = u(i), where 
v(t) is the solution with initial value v(0) = vq. Using (4.32) we obtain the following 
corollary 

Corollary 5.4. TgSnF(t) = SoF(t)Tg. 

It follows that periodicity of trajectories is conserved by (4.9). 

Proposition 5.5. If v(t) is the solution of (4.9) and v(0) is a periodic trajectory 
with period T then v(t) is periodic with period T for any t. In particular, ifv(0) is 
constant trajectory, then v(t) is constant for any t. 

Recall that for v(t) £ X we denote by v(t, s) the value of v(t) at the parameter 

sel. 

Lemma 5.6. Suppose v is a solution of (4.9) such that \\v(0)\\x < 2vkqG and 

(5.9) v(t, s) = v(0,t + s) for all i.sgR. 
Then u(t) = v(0,t) + W(t, v(0, •)) is the solution of (2.5) satisfying 

P n u(t) = v(0, t) = v(t, 0) for all tel. 
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Proof. We first note that by (5.9) and (4.32) 

W(0; v(t, •)) = W(0; T t v(0, •)) = T t W(0; v(0, •)) = W(t; v(0, ■)) , 
so we can set 

(5.10) w(t) = W(0; v(t)) = W(t; v(0)) . 

Since ||u(0)|| < 2vkqG we have <p = l, where ip is as in (4.10). By (5.9), (5.10), and 
the definition of the map W, we have 



(5.11) 



as well as 



(5.12) 



^-w(t) + vAw{t) + Q N B(v(0, t) + w(t)) = 
at 

^-w(t) + vAw{t) + Q N B(v(t, 0) + W(t; v(0, •))) = Q N f 
at 



j t v{t) + vAv{t) + P N B(v(t) + w(t)) = 
jv{t, 0) + uAv(t, 0) + P N B(v(t, 0) + W(0; v(t, ■))) = P N f 



Adding the equations (5.12) and (5.11), we see that u(t) = v(0,t) + W(t,v(0, •)) is 
a solution of (2.5). □ 

We also have the following property of (4.9) 

Proposition 5.7. Suppose that Ui(0, s) = V2{0,s) for s < sq. Then vi(t,s) = 
V2(t, s) for s < so for any t. 

Proof. By Corollary 5.1 the map W is well-defined on the space C&((— oo, so] , PnH). 
Thus we can consider the differential equation (4.9) on this space. It is still 
an ordinary differential equation, and so through each point there is only one 
solution. Each solution of (4.9) on the whole X gives a solution of (4.9) on 
Ci,((— oo, so], PnH), by restriction. Since solutions vi(t) and V2(t) of the restric- 
tion of (4.9) on Cb((—oo, sq], PnH) coincide for t = 0, they coincide for all t. □ 

6. Absorbing ball 

We will now prove that the determining form (4.9) has an absorbing ball in X. 
Recall that this means an estimate of the type (2.24), holds but with || • || replaced 
by | ■ \x- Let g = Pxf and take the scalar product with v = v(t,s) (and denote 
W(v) by w) to obtain 

\j t \ v ? + HMI 2 < IflIM + \<p{B{v,v),w)\ + \ip(B(w,v),w)\ 

< \g\\\v\\ + CBfWvl^lwlihieN) 1 ^ 2 + ci,tp\w\\\w\\\\v\\ 

<r abU II . II »2 GP(G,N) n 1/a 2 G 2 f3(G, iV) 2 

< hqis — G\\v\\ + vcb\\v\\ — (IneA) ' + k q v c l — ||u||, 

where we used inequalities (2. 31), (2. 33), (4.24) and (4.28). It follows from Lemma 3.1 
that if 

(6.1) N > c B G/3(G, N)(\neN) 1/2 

holds, we have the following 
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Proposition 6.1. If (6.1) and (5.6) hold, then 
(6.2) Id < vG 



Is] | rrr l3(G,NY 



While the relation ||i>|| < Ko-^M provides an immediate bound on ||w|| provided 
(6.1) holds, we can derive a sharper bound on \\v\\. As above, taking the scalar 
product with Av we have 
(6.3) 

--|M| 2 + HAi>| 2 = (g,Av) - (B(v + w,v + w),Av) 

< \g\\Av\ - (B(v,w),Av) - (B(w,v),Av) - (B(w,w),Av) 



< \g\\Av\+c T (2\\v\\\\w\\\Av\ + ||^|| 2 |^|)(ln eA^) 1 / 2 

j^G + c T (lne7V) 1 / 2 G 2 (6 + /?(G, N))(3{G, N) 



^22 



\Av\ 



where v = v(t, s) and (6.3) is valid at each set. Now applying Lemma 3.1 (i) and 
then (ii), we obtain 

Theorem 6.2. // (5.6) holds, then the ball centered at the origin, with radius 
(strictly) greater than 



+ c T (lneiV) 1 / 2 G(6 + /3(G, N))0(G, N) 



(6.4) r = vkqG 
is an absorbing ball in X for (4.9). 

Remark 6.1. The 0(G 2 ) estimate (up to a logarithm) for the radius of the absorb- 
ing ball for (4.9) in (6.4) is significantly larger than the 0{G) estimate for radius of 
the absorbing ball for the NSE. These estimates can be brought closer if we assume 
N is larger. Indeed, if we take for some 7 > 1 

(6.5) N > (6c T G(lneiV) 1/2 ) 7 , 
then 

mN) < 9^/(6^) + 1 < 3NU2 < ^ 
N-N~/2 N 

and consequently 

G(lneAT) 1 / 2 ^(G,A^) < — N^r . 

6ct 

Thus, if equality holds in (6.5) with 7 = 3/2, we have (up to a logarithm) that 
?'o = 0(G 3 / 2 ), and if (6.5) holds with 7 = 2, we have ro = 0(G). 

Let B r be an absorbing ball in X, centered at with radius r > ro, and let 
•M = r\t>oSr)F{t)B r . Note that M is independent of r: vq € X belongs to M if and 
only if there is a bounded solution v : (—00, 00) — > X of (5.3) defined on the whole 
real line such that v(0) = vq. We can also characterize M. as the maximal bounded 
set that is invariant (i.e. S^(t)J^A = A4 Vt > 0), both forward and backward in 
time, under the flow generated by (4.9). We do not know if M. is the analog for the 
determining form (4.9) of the global attractor of the Navier-Stokcs equations (see 
discussion leading to (2.28)). 

By Theorem 6.2 M is a subset of B ro . 
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Lemma 6.3. If v is a bounded solution of (5.3) defined on the whole real line and 
\\W(v(t))\\Y < vkq y/^G/cL for all t then \\v(t)\\x < SvkqG for all t. In particular 
this is true ifN> V2cI(9c T G 3/2 (ln eN) 1 ' 2 + G 1 / 2 ) and N > 6c L G. 

Proof. Let w(t) = W(v(t)) and w(t,s) = w(t)(s). Suppose that \\v(t', s')\\ > 3vk G 
for some t' and s' . Then since v is continuous there arc to < t' < t\ suth that in the 
interval [ioi^i] we have \\v(t, s')\\ > SvkqG. On this interval the map t <— > v(t,s') 
satisfies the equation 

j t v(t, s') + vAv{t, s') + B(w(t, s')) = g. 

Taking its scalar product with Av(t, s') we obtain inequality 

j t \\v(t, s')\\ 2 + v\Av(t, s')\ 2 < \g\\Av(t, s')\ + \(B(w(t, s'), Av(t, s% w(t, s'))\ 

and using (2.33) we have that \(B(w, Av),w)\ < c h \w\\\ w\\ \A 3 / 2 v\ < c h \\w\\ Y \Av\, 
and hence 

j f \\v(t, s')\\ 2 + v\Av(t, S ')| 2 < (<Gv 2 k 2 + c h \\w\\ 2 Y )\Av(t, s')\. 

Note that by Lemma 3.1 (i) to = — oo, since Giskq + c-L\\w(t)\\ 2 / (vko) < ZGvkq, but 
by the Lemma 3.1 (iii) this is impossible. □ 

Remark 6.2. It is easy to give an analytic characterization of the solutions of (4.9) 
in M. . Indeed we have 



v{t) = e-^-^vito) + / e~ u{t - e)A [g - P N B{v(9) + W(v(8)))]d0, 
and letting to — > — oo we easily obtain the relation 
(6.6) v{t)= f e~ v(t - B)A [g-P N B{v-{e) + W{v{e)))]d8 t > 0, 



in X. Clearly (6.6) implies that v(t) is in M. 
In the case when 

v(t,s) = P N u(t + s), s,teR, 

where u(t) is a solution of (2.5) in A, it is easy to see that the relation (6.6) is 
equivalent to 

(6.7) «(*)=/" e-^ s - a)A [g-P N B(v(a) + W(v(a)))]da set. 

J — oo 

It follows that 

v(s) = Pnu(s), set, 

where u(s) is a solution of (2.5) in A if and only ifvGX and satisfies (in H) the 
relation (6.7). 
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7. Traveling waves 



We say that v(t, s) = SBF(t)vo(s) is a traveling wave solution for (4.9) if 



Denote w(t,s) = W(s;v(t, •)), u(t,s) = v(t,s) + w(t,s), and define iuq € Y, u <G 
CbQR, H) by wo(s) = w(0, s), ito(s) = u(0, s). Suppose u is a traveling wave solution. 
Then by (7.1), (4.32) w(t, s) = wo(t + s), thus u(t, s) = uo{t + s). We will provide an 
equation for uq and use it to derive the properties of translation invariant solutions. 
Using (7.1), equation (4.9) can be written as 

(7.2) — vo + vAv = P N f -P N B(<pv +w ), 
at 

while (4.13) is written as 

(7.3) ^-w + vAwq = Q N f - Q N B(tpv + w ). 
as 

Changing t to s in (7.2) and combining the result with (7.3), we have the equation 

for Uq 

(7.4) — u + pAu = f -B(<pv + w ). 
as 

Taking the scalar product with Auq, we have 



Using (2.15) and (2.13), we obtain 
{B(ipv + w ), Au ) = ip[(B(v , w ), Av ) + (B(w ,v Q ), Av )] + {B(w , w ), Av ) 



(7.1) 



v(t, s) = v (t + s). 



(7.5) 



^ll«oll a + H^«ol 



(f,Au ) - {B(ipv + w ),Au )- 



i 2 (B(vo,v ),Aw ) + (p[(B(v ,w ),Awo) + (B(w ,v ),Awo)] 
= (V 2 ~ l p)(B(v ,vo),Aw ) + (1 - (p)(B(w ,w ),Avo) , 



so by applying (2.15) again, we have 

\{B(<pv +wo),Au )\ < (p\(B(v ,v ),Aw )\ + \(B{w ,w ), Av„)\ 

< ip\(B(v ,w ),Av )\ + <p\(B(w ,Vq),Av )\ + \(B(w ,w ), Av a )\ . 
We now apply (2.30), (4.16), (4.31) to obtain 



--\\u \\ 2 + v\Au \ 2 < \f\\Auo\+c T (Q + p{G,N))v 2 KlG 2 p{G, N){\neNff 2 \Avo\ 
< v 2 k 2 G[1 + c T {6 + f3(G, N))G(3(G, AT) (In eN)^ 2 ]\Au \ 



and consequently by Lemma 3.1 (iii) 

(7.6) ||uo(s)|| < vk G[1 + c T (6 + 0(G, N))G(3(G, N)(ln eN) 1 / 2 ] . 
We next calculate an absolute constant c' T such that if 

(7.7) N > c' T Gj3(G, N) In eN , 



then 



c T (6 + /3(G,^))G/3(G,7V)(lneiV) 1 / 2 < 1 , 




IKOOII < ll«o(s)ll <2vk G. 
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First note that for c' T > 18ct, we have /3(G, N) < 2, for all N > 1. Then estimate 

» rMr Ai\n ArtVa s Q^GlneiV + (lneiV) 1 / 2 ^!/!/! 
8c T G(3(G,N)(lneN) = 8c T ^ _ ejv)1/2 

9c T N/c' T + N/cL 72c 2 T + 8c T 
N-3c T N/c' T c' T -3c T ' 

so that it suffices in (7.7) to take 

(7.9) c' T > 72c| + IIct ■ 

Theorem 7.1. Suppose (7.7) holds for c' T as in (7.9). Then a solution v(t) of (4.9) 
is traveling wave solution if and only if v(0) is a projection of a solution of (2.5) 
in A. 

Proof. The fact that any traveling wave solution of (4.9) is the projection of a 
solution of (2.5) follows from (7.8) and Lemma 5.6. To prove the converse, recall 
that any solution u(t) of the (2.5) which lies on the attractor can be extended to 
a solution to the complexified NSE which is analytic and bounded on some strip 
< c function into H (see e.g. [4]). Let u(t) be a solution of (2.5) defined and 
bounded on the whole real line. Consider a curve in X, t n- (s n- p(t + s)), where 
p = Pnu. We need to show that this curve is differentiable and its derivative is 
s 1— > dp/dt(t + s); thus we need to prove that for any real t 

(7.10) sup|| P(* + At + .)- P (t + 5 ) _ d (f + a)|| _ as A ^ . 

s6R At ">t 

Since p is analytic and bounded in the strip < c, its second derivative is also 
bounded in a smaller strip. The limit in (7.10) now follows from Taylor's theorem, 
specifically 

n p(t + At + s) - p(t + s) d . ... . „ d 2 , 

sup|R ^ i-_p(t + s )|| <su P Ai sup II— p(r + «)|| . 



seR Ai sS r re[t,t+t+A*] 



□ 



Remark 7.1. Note that a function v(t,s ~ t) being independent of t is equivalent 
to v{t,s) = f(0, s + t), for all s,t; thus if \\v(0)\\x < 2vkqG, then v(t) is the Pn- 
projection of a solution on the global attractor of the Navier-Stokes equations. 

8. Equicontinuity of W 
We have from (4.17) and (4.26) that 

f \Aw\ 2 < E^S! - s ) + E , 

for some constants Eq,E\. We also have for any u E H 

\(B(ipv + w), u)\ < 9c B G 2 u 2 K 2 Q (\neN)^ 2 \u\ + 3(c A + c T )v 2 n 2 G 2 /3(G, JV)(ln eA^) 1/2 |u| 

+ c A v Ko ^[3(G,N)\Aw\\u\ ; 

thus 

(8.2) \B(<pv + vj)\ < E 2 +E s \Aw\ 



(8.1) 
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where 

E 2 ={9c B + 3(c B + c T )j3(G, N)]v 2 K 2 G 2 p{G, A^)(lneiV) 1/2 
Integrating (4.13) from so to si we have 

/•si rsi nai 

\w(si) - w(s )\ < v j \Aw\ + I \B((pv + w,(pv + w)\ + / \h\ 

J s J s J s 

<(v + E 3 )( Sl - s ) 1/2 [ J Bi(s 1 - so) + E } 1/2 + (E 2 + \h\){ Sl - s ) 

We have proved the following 

Lemma 8.1. The family {W{v)(-)\v <G Cf,(R, PnH)} is equicontinuous on R. 

Let v(t) be a solution of (4.9), and denote w(s,t) = W(v(t))(s). Then by 
Lemma 5.2 



\\w(s,t) - w(s,t )\\ < L w \\v{t) - v{t )\ 



Since v(t) is a solution of an ODE, it is continuous. Thus for any finite interval 
[— T, T] there is a constant Ct, such that 

||«(t)-v(fo)||x<Cr|t-to|, t,t €[-T,T}. 

This proves 

Proposition 8.2. For any solution v(t) of the equation (4.9) the function w(s,t) = 
W(v(t))(s) is equicontinuous on R 2 . Moreover, there are constants c, c\ and Ct 
such that 

\w(s,t) - w(s ,t Q )\ < c\s - s | 1/2 + a\s - s | + —\t - t \, 

whenever t, to £ [— T, T] . 

We now assume that (5.6) holds, and that h G D(A 1 ^ 2 ). Then by (4.36) we have 
(8.3) \AW(v)(s)\ < T; , for all v e X, s e R , 

where T is the right handside of (4.36). We have therefore that 

\\w(s,t)-w(s ,t )\\ < (2T)^ 2 \w(s,t)~w(s ,to)\ 1/2 



(8.4) r ^_ iV2 

< (2L) 1 / 2 



c 

c\s - s \ 1/2 +d\a- s | + -£\t — 1 | 



iV 

for all s , si e R, and all i, t e [-T, T]. 

Recall the extension of the Arzela-Ascoli theorem to the case of Banach-valucd 
functions. 

Theorem 8.3. Let K be a compact metric space, X a Banach space, $„ : k t— > 
<&n(k) € A", iwif/i {$n}^Li equicontinuous and {& n (k)}^ =1 relatively compact for 
each k 6 A'. TTien t/iere exists a subsequence {& nj }j C Li which converges in C{K; X). 



22 



C. FOIAS, M.S. JOLLY, R. KRAVCHENKO, AND E. S. TITI 



9. Weak limit of time translated solutions 

According to Remark 7.1 (Section 7), a solution v to the determining form (4.9) 
in the set M. is the PAr-projection of a solution of the Navier-Stokes equations (2.5) 
in its global attractor A if and only if 

(9.1) \\v(t, s)\\ < 2vkIG , Vt,seR, 
and 

v(t + t n> s-(t + t n )) = v(0, s) , Vf.sel, 

for any sequence {t n }. It is therefore instructive to study the asymptotic behavior 
of a sequence 

(9.2) v n (t,s)=v(t + t n ,s-(t + t n )) , n = l,2,... 

where t n — s- oo, and v(-) £ M. is a solution of (4.9) satisfying (9.1). We have the 
following lemma, 

Lemma 9.1. There is a subsequence {rij} and functions v, id and R such that 
v,R£ C(R, L°°(R, H) weak *), w £ C(R 2 , V) and the following hold for each T, S 

(9.3) Urn sup IK (*,•)-«(*,•), 01=0 for^£L\R,H) 

i te[-T,T] 

(9.4) lim sup \(B(v nj (t,-)-v(t,-))+R(t,-),0\ =0 for^£L\R,H) 

i te[-T,T] 

(9.5) w nj (t,s) = W(T tn .v)(s-(t + t nj ))^w(t,s) (Mel) 
in C([-T,T] x [S,S];V), for all T,S>0. 

Proof. It follows from (4.9) and (8.3) that there exists a constant C 

dv n (t,s) 

I — Q- f — \<C (yt,s,n) 

Choose a dense subset in the unit ball of L 1 (R, H). Then it is well-known 
and easy to check that 

^ 2 fc (l + \{p-p',Ck)\) 

defines a distance in the ball of any fixed size in L°°(M.,H). Let us compute 
d(v n (t,-),v n (t',-)). Note first that 



(v n (t,s) - v n (t',s),£ k (s)) ds\ 

= | J (v(t + t n , S - t - t„) - V(t' +t n ,S - t' - t n ), ^fc(s)) ds\ 

<C|t-t'| + l / (v(t + t n ,s-t-t n )-v(t + t n ,s-t' -t n ),£ k (s)) ds\ 



= C\t-t'\ + | / (v(t + t n ,s-t n ),Zk(s + t) -Z k (s + t'))\ ds 

< C\t-t'\+3GvKoU k -n-t>£k\\v 

which, provided k is fixed, goes to when t — t' goes to 0, uniformly in n. It follows 
that for any fixed e there is a 5 such that \t— 1'\ < 5 implies that d(v n (t, •), v n (t', •)) < 
e for all n. This means that the set of functions t \— > v n (t, •) is an equicontinuous 
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family of functions in C([—T, T],L°°(R, -ff) wea k*) for any T. Note that \v n {t,s)\ < 
\\v\\x/ko, therefore v n (t,-) belongs to the ball of radius ||u||x/ko hi L°°(R 7 H), 
which by the Banach-Alaoglu Theorem is compact in L°°(R, iJ) wca k* . Using the 
Arzela-Ascoli Theorem 8.3 and the diagonal process, we obtain a function v in 
C(R, L°°(R, i?)weak*) such that for some subsequence of {v nj }, (9.3) holds for any 
T. Without loss of generality we assume that (9.3) holds for the sequence {v n }. 

To deal with the sequence R n {t, •) = —B(v n (t, •) — v(t, •)) note that it is a com- 
position of an equicontinuous family {v n (t, •) — v(t, •)} in C(R, L°°(R, -P/v-ff ) WG ak* ) 
and a map L°°(R, P/v H) wea k* — > L°°(R, 7J) we ak* induced by B. Since PnH is finite 
dimensional, B is a continuous quadratic form on PnH, and hence the induced 
map is continuous. It follows that the family R n (t, •) <G C(R, L°°(R, -ff) wea k*) is also 
equicontinuous. Now note that by (2.35) 

\B{v n (t,s)-v(t,s))\ < c A \v n (t,s) - v(s)\ 1/2 \Av n (t,s) - Av(s)\ 1/2 \\v n (t, s) - v(s)\\ 
< c A (K N) 2 \v n (t,s)-v(s)\ 2 < 2c A ( Ko N) 2 \\v\\ x , 

thus R„ (t , • ) is in the ball of radius 2ca (ko N) 2 \ \ v \ \ x in L°° (R, H ) . Therefore, again 
by Theorem 8.3 and the diagonal procedure, there exists R(t, •) £ C(R, L°°(R, if) wca k* 
and a subsequence such that (9.4) holds for any T. Once again, without loss of gen- 
erality, we assume that (9.3) and (9.4) hold for the whole sequence. 

Finally, w n is equicontinuous in C(R 2 ,V) by Proposition 8.2, and for each t,s 
w n (t,s) is precompact in V by (8.3), thus, by Theorem 8.3, (9.5) holds for some 
subsequence v nj . □ 

The main result in this section is the following. 

Proposition 9.2. Suppose that N and G satisfy (5.6). Let v, w and R be as 
in Lemma 9.1. If v(t,s) and w(t,s) are both independent of t, then R is also 
independent of t and 

u(s) = v(s) + w(s) 
satisfies the Navier-Stokes equation 

(9.6) -r-w(s) + vAu(s) +B(u(s),u(s))) = f + R(s) , seR. 

as 

The function R can be viewed as the analog of the Reynolds stress due to aver- 
aging the Navier-Stokes equations in the theory of turbulence [15, 13] Thus, u(-) is 
a solution to the Navier-Stokes equations (2.5) with the original force /, if and only 
if v(t, s) and w(t, s) are independent of t, and R(s) = 0. We do not know if any of 
these conditions are superfluous. 

Proof. Using a change of variables, we write 
(9.7) 

w n (t,s) = w(t + t n ,s- (t + t n )) = W(r tn v)(s - (t + t n )) 

fS-(t-t„) 

e - v{ s-t-t n -a)A ^ _ Q N B(v(t + t n , a) + w(t + t n , a))} da 
{vA^h- f e-^ s -^ A {Q N B(v n (t,a)+w n (t,a))} da. 
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Now write 
(9.8) 

w(s) - (uA)~ 1 h + / e- u( - s -^ A Q N [B(v(a)+w(a))-R{t,(r)]da 



w(s) - w n (t, s)- 

s 

e -u(s-a)A QN [Bfafa CT ) + Wn ( t: _ B (v{ a ) + w(a)) - B(v n (t, a) - v(a))} da 
Notice that 

B(v n + w n ) - B(v + w) - B(v n - v) = 

(9.9) 

B sym (v n + w, w n - w) + B(w n -w) + B sym (v + w,v n - v) 

where 

B sym (u 1 ,u 2 ) = B(ui,u 2 ) + B(u 2 , mi) . 
Using (2.35), (8.3) and (2.34) as well, we find that 

\B(v n +w,w n - w)\ < c A \v n + w\ 1/2 \A(v n + w)\ 1/2 \\w n - w\\ 

<—(\Av n \ + \Aw\)\\w n -w\\ 

< C A (n\\v\\ X + £j ||«>n-tS|| . 



\B{ Wn -W,V n +w)\< C L \w n - W\ 1/2 \\w n - H| 1/2 K + w\\ 1/2 \MVn + w)\ 1/2 

<^^(\Av n \ 2 + \Aw\ 2 ) 1/2 \\ Wn -w\\ 
- ~ N i/2 ( K o N \\ v n\\x + \Aw\) \\w n - w\\ 
<^T72( K oN\\v\\ x + r)\\w n -w\\ , 

\B(w n - w)\ < c A \w n - w\ 1/2 \A(w n - w)\ 1/2 \\w n - w\\ 

< —^-\A(w n - w)\\\w n - w\\ 

K N 

< —% T W W n - W\\ , 
K N 

so that for each t 

- v ( s ~ a ) A Q N ]J3 syra (y n + w , w n — w) + B(w n — w)] da — > as n — > oo , 



e 

— oo 



by the Lebesgue dominated convergence theorem. Indeed, the integrand is bounded 

by 

(9.10) 2f|H| y e-^ s - CT > K « 
where 

f = ( CiW + c t ^)|HU + (^ + -J I7i + ^|)r, 

and goes to at each a since ||u>„(i, a) — u>(cr)|| — > when n — > oo. 
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Let k € D(A). Note that the function a ^ e~ v( - s -^ A k if a < s and a i-> if 
cr > s is in L 1 (M, iJ). Thus by (9.4) 



(9.11) / e" l/(;i -' T)A (Qjv(B(u„(t,cr)-v(cr))+i?(t,cr)),fc) dcr^O as n ^ oo . 

J — oo 

We also have 

(/ e-^ s -^ A Q N B(v + iB,v n -v) da,k) 

Q N (B(v + w, e- u{s - a)A k),v n -v)da. 

But by (2.35) we have that, as a function of a 

B(v(a) +w(a),e" j{s - a)A k) e ^((-oo, s], H) . 

We conclude that 

^ e -u( a -a)AQ NB ^j + ^ Vn - y) 7 k) da ^ as n -> oo . 

For the other half of B sym (v + w, v n — v) wc write 

(e-^ s -^ A B(v n -v,v + w), k) = (v n - v, 0(a)) 
where (due to (2.36)) we have 

\6{a)\ < \\v(a) + w(a)\\e-^ a -^ Ko \k\ 1/2 \Ak\ 1/2 V a < s . 
It follows that 0(a) € L°°((— oo, s]), and hence 

(e~ u( - s ~ a)A Q N B(v n - v, v + w), k) da -> as n -> oo . 

In particular, for fc € D (A), (9.8) now implies that 

(w(s) - (iM)"^ + f e-^ s -^ A Q N (B(v(a) + w(ct)) - R(t, a))da, k) = 

J — oo 

Since -D(A) is dense in iJ, we deduce that for each t 

(9.12) w(s) = (vA)- 1 h- [ e-^ s -^ A Q N (B(v(a)+w(a))^R(t,a))da . 



Note that (9.12) implies that QjyR(t, s) does not depend on t. Indeed, by the vector- 
valued version of the Lebesguc differentiation theorem we obtain that for each t 
QNR(t,s) = r(s) for almost all s, fo some function r G C(M, L°°(R, Qjv-ff)weak*) 5 
which means that Qat-R is constant as a function in C(R, L°°(R, -ff) wea k* )■ 
Now observe that 

v n (t, s) = v(t + t n ,s -t- t n )) 

rt+tn 

e -u{t+t n -r)A^ g _ p N B( v (r, s-t-t n ) + w(t, s-t- t n ))} dT 
e -v(t-r)A^ g _ PnB ( v ( t + tn ,s-t- t n ) + w(t + t n , s-t- tn))] dr 

»(t-T)Afe _ P N B( Vn (T, S + T-t)+ W n (T, S + T- t))} dT 



— OO 

— oo 
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for all i,seK; consequently 
(9.13) 



S-2 



(v n (t,s),k) ds = ((uA)- 1 g,k)(S 2 -S 1 ) 
>s 2 



{B(v n (T, 8 + T-t) + W n (T, S + T- *)), e^^P^fe) ds 



dr 



=((vA)- 1 5,fc)(S , 2 -5i) 

»S 2 +T-t 



(B(w„(t, s) + u; n (r, «)), e-^-^ A P N k) ds 



Si+r-f 



where fc € -D(j4), Si < S2 are fixed (but otherwise arbitrary). By the convergence 
in (9.4) we can now pass to the limit in (9.13) by going through steps analogous to 
those from (9.7) to (9.12). Indeed, if k is an eigenvector of A with Ak = Xk, then 
by (9.9) 



r S 2 +T-t 


•/si4 


(B(v n (r, s) + W n 




T-t 




rS 2 +T-t 


r 


/ (B(v + w)- 


—00 


JSl+T-t 




[■S 2 +T-t 


r 


I {B svm {v n + 


—00 


JSl+T-t 



dr 



R,e-" (t -^ A P N k) ds 



dr 



dr 



-v(t—r)\ 



o -v(t- T )X 



S-2+T-t 

Sl+T-t 
S 2 +T-t 

Sl+T-t 



(B sym (v + w,v n - v),P N k) ds 



dr 



(R + B(v n - v),P N k) ds 



eh 



On the right-hand side the second, third and fourth terms go to as n — > 00. 
In the second term we use a bound similar to (9.10) and the Lebesgue dominated 
convergence theorem. For the third and fourth note first that the integrands are 
bounded by L 1 functions, thus it suffices to show that for fixed r the integrals 



S 2 +T-t 



Si+r-i 



(B sym (v + w, v n - v), k) ds 



and 



S 2 +T-t 



St+T-t 



(R + B(v„ -v),k) ds 



go to as n — > 00. For the latter, note that the function X[s 2 +T-t,s 1 +T-t](s)k is in 
L 1 (K, H), where xi is the characteristic function of an interval 7, and use (9.4). For 
the former, write (B sym (v + w,v n ~ v),k) — (v n (t,s) — v(s),ip(s)) and note that, 
as above, by (2.35) and (2.36) X[S 2 +T-t.s 1 +T-t]( s )4'( s ) is m i 1 (R, H) and therefore 
we obtain the convergence by (9.3). 
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Setting u = v + w, we obtain 

r-S 2 

(v(s),k) da-dvAy^k^Si-S!) 

»S 2 +T-t 

Sl+T-t 
S 2 



Si 



Si 



Si 



(B(u(s)) - R{T,s),e- u{t ~ T)A P N k) ds 
{B(u(s +T-t))- R{t, s + t- t),e- u{t - T ^> A P N k) ds 
e- u{t - T)A P N [B{u{s +r-t))- R(t, s + t- t)] dr, k) ds 
e-< s -^ A P N [B{u{a)) -R(t-s + a, a)] daX ds, k) . 



dr 



Since for fixed Si < S2, t in R, k can run over all linear combinations of eigen- 
vectors of A, of which arc dense in H, we deduce 



S-2 



v(s) 



-v{s—<j)A 



[g - P N B{u{a)) + P N R(t - s + a, a)} da \ ds = 



for all Si < S 2 , t in R. 

The final function in braces above is in S, S], H) for all S > 0. Therefore, 

the vector-valued version of Lebcsgue's differentiation theorem now implies 



v(s) = ( e- v{s - a)A [g - P N B{u{a)) + P N R(t - s + a, a) dt] da 

J —OO 



(9.14) 



for s £ Et, where the measure of E t is 0. Note that since C(R, L°°(R, i/)) W cak* C 
L 1 1 oc (IR 2 ), both the right- and left-hand sides of (9.14) are in i, 1 oc (]R 2 ), and therefore 
by Fubini's theorem (9.14) holds for almost all pairs (t,s) £ R 2 . By substituting 
(t, s) with (t + r, s + t) we obtain that for almost all t, s 



(9.15) v(s + t) 



o -v(s + T-(j)A 



[g - P N B(u(a)) + P N R{t - s + a, a) dt] da 



holds for almost all r. Fix such t, s that (9.15) holds for almost all r. Then the 
right-hand side of (9.15) is continuous in r and hence v coincides almost everywhere 
with continuous function. Thus, we can suppose that v is (absolutely) continuous. 
Taking the derivative with respect to r of (9.15), we obtain that for almost all 

{t,s) £ M 2 , 



dv(s 



ds 



+ vAv(s + t) + P n B(u(s + t),u(s + t)) = g + P N R(t + T,s + t) , 



for almost all r. It follows that 

+ vAv{s) + P N B{u{s),u{s)) = g + P N R{t, s) , 

ds 

for almost all t, s. We therefore obtain that the function R £ C(R, L°°(R, i?) W eak*) 
is almost everywhere equal to some fixed element of L°°(M.,H). Since R is contin- 
uous, this equality holds everywhere, thus R docs not depend on t. □ 
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10. Stationary solutions of (4.9) 
Recall that uq £ H is stationary for (2.5) if 

vAu + B(u ,u ) = f. 

Clearly, if 

(10.1) w(O,0 =PnM-) > 

then v(t, •) — PnUo(-) for all t is a stationary solution of (4.9). In fact by Theorem 
7.1 any stationary solution v(t, s) of (4.9) which is constant in s must satisfy (10.1). 
The set of all stationary solutions of (4.9) may, however, be much larger than the 
set of constant functions in X that are the projections of the stationary solutions 
of (2.5) as in (10.1). In fact we show in [11] that a determining form for the Lorenz 
system has stationary solutions which are not in the global attractor of the Lorenz 
system. 

A solution v(t, •) is a stationary solution of (4.9) if and only if v(t, s) = vo(s) for 
all t and 

uAv {s) = g- P N B{ V v a (s) + W(v )(s)), 
that is, if and only if it satisfies the following differential algebraic equation 

(10.2) vAv (s) =g- P N B(pv (s) + w (s)), 

(10.3) + uAw (s) =h- Q N B(^v (s) + w {s)) , 
as 

where wq = W(vq). In this case the energy and enstrophy balances are 

(10.4) ~\w \ 2 + u\\u \\ 2 = (f )Uo ), 

(10.5) ~\\w Q \\ 2 + is\Au Q \ 2 = (f,Au ), 

where uq = vq + wo (compare to (2.20), (2.21)). 

It is shown in [8] that solutions on the global attractor of the NSE must satisfy 
HKOII 2 < l/IN*)l for a11 t e R, with equality holding for some t if and only if / 
is an eigenvector of A. One might ask if at least the stationary solutions of (4.9) 
satisfy an analogous property. We will give below a partial answer for this case, 
which involves an application of (10.4), (10.5). 

Let wo = vq + wq . We first prove the following: 

Lemma 10.1. Let wq and vq satisfy (10.2) and (10.3), and let oo < a\. Suppose 
that |wo(s)| is nondecreasing in (ao,ai), then 

(i) ^||u (.s)|| 2 < |/||u (s)| for all s G (<r ,ai). 

(ii) If, in addition, u>o{s) is nonzero for all s € (<7o,<7i), ||u;o(s)|| 2 /|u'o(s)| is 
nondecreasing, and for some so G (oo,o~i), we have ||tto(so)|| 2 /|wo(so)| = 
|/|/V, then f is an eigenvector of A for the eigenvalue A = |jwo(so)|| 2 /|uo(so)h 
A > (A^o) 2 , andf = Q N f. 

Proof. In the proof we drop the subscript on wq, vq, uq, to simplify notation. Denote 
Xw = \\w\\ 2 /\w\, X u = \\u\\ 2 /\u\, X w = \\w\\ 2 /\w\ 2 , X u = \\u\\ 2 /\u\ 2 . Since w(s) ? 
in (<Jo,cri), these quantities are well-defined. From (10.4) and the fact that \w\ 
is nondecreasing it follows that (f,u) — v\\u\\ 2 > 0, which in turn implies that 
Xu(s) < \f\/f, proving (i). 
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Passing to the proof of (ii), we start by observing that from (10.4), (10.5) it 
follows that 



Xw{s) ~ Xw(s ) 



(10.6) 



2v 

N 



\i\ 2 -4- 



'2v 



A u > — A u ( : f 
v 



da. 



By (i), since \w\ is nondecreasing, and the assumption Xu(so) = l/l/^i we deduce 
from (10.4) that (/, u) — \ f\\u\\ hence we have u(s Q ) = cf, where c = |/| 2 /H|/|| 2 - 

Since Xw is assumed to be nondecreasing, the integrand in (10.6) is nonnegative 
for almost all a; also every summand in it is continuous, except for |//2// — (A — 
A u /2)u|, which is lower semicontinuous. This suffices however, to ensure that the 
integrand is in fact nonnegative at every point of (co, 0*1 )• 

In particular, at the point so we have that ||m|| 2 = (f /v, u), thus it follows that 

A,, 



(10.7) 



is- 



■*(so) 



u(s )\ = 0, 



hence (A - X°)f = 0, where A 



A u ( So )- Indeed, 



If A < (Nkq) 2 it would follow that / = g, which gives that w(sq) — 0, but this is 
contrary to our assumption. Thus we have that A > (Nkq) 2 , and so / = h, g = 0, 
v{s ) = 0, w(s Q ) = h/v\° and Ah = X°h. □ 

Let h = h/v\ Q . Now we need the following lemma: 

Lemma 10.2. Assume that vq and wq satisfy (10.2) and (10.3), and that g = 
and h is an eigenvector of A. Then for large enough N, v(s) — and w(s) = ho 
for all s, where A° is an eigenvalue of A. 

Proof. Once again, denote Wq,v 0i uo by w, v, u for simplicity. From (10.2), (2.32), 
and (4.28) we obtain that 

HM| 2 = ~(B{u,w),v) = —(B(w,w),v) + {B(v,v),w) 



whence 
(10.8) 



v\\v\\ 



< (B(w,w),v) +ct\\v\\ 2 \\w\ 



ct ; m eN ) ' 



(lneAO 1 ^ 



k N 



< 



(B(w,w),v). 



(10.9) 



By denoting w = h a + S we have 

-(B(w, w),v)=- (B(S, w), v) - (B(h Q ,S), v) 

< (In eN)^ 2 \\S\\(c T \w\\\v\\ + c B \h \\\v\\) 



< 



1/a vGP(G,N) 
(IneiV) ' (c T — 



c B v 



^)ll«llll«ll 



From (10.8) and (10.9) we obtain 
(10.10) 

Gj3(G, N) 



1 — ct- 



N 



-(IneiV) 1 / 2 



< (ln e 7V) 1 /2^( CT/3(G ,iV) + Cj 



kAN, 



A° 
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Therefore if 

(10.11) N > 2c T G(3(G, N)(\n eN) 1 / 2 , 
we have the bound 

(10.12) \\v\\ < ci\\8\\, 
where 

Cl = 2(lneN)^^(c T t3(G,N)+c B ±). 
Now using (10.3) and (10.12) we obtain that 

~\S\ 2 + v\\6\\ 2 = -(B(u,u),5) = -(B(v + h + S,v + h ),S) 
I as 

= -(B(v + S, h ),6)-(B(u,v),S) 

<Qr(]neA0 1/2 |MIIMII*l + c L |*|||<f|||NI + crQneAO^INIMII*! 

< v{c TCl {lneN)^ + ^\\5f + ^(IneiV)^ Mx + ^{G,N) 2 _ 

N A kqN 

It follows from Lemma 3.1 and Theorem 6.2, that if N is such that (10.11) holds 
and moreover 

N 2 > 2G(c T ci(lneiV) 1 / 2 + c L ), 

(10.13) 

TV > 2c TCl (G + c T (lneA0 1/2 G 2 (6 + /3(G, N))/3(G, N) + Gf3(G, N)), 

then S(s) = for all s. Therefore from (10.12) we have for all s that v(s) = , and 
hence u(s) = h/v\° = w(s). Since h is an eigenvector of A, it follows from (10.4) 
that Ah = X°h □ 

Therefore we obtain the following proposition. 

Proposition 10.3. Suppose N satisfies (10.11) and (10.13). If the assumptions 
of Lemma 10.1 hold for wq = W^fo) for some stationary solution vq of (4.9), then 
any stationary solution of (4.9) is zero. 

Wc finish by simplifying the conditions on N in Proposition 10.3. Suppose that 
G > 1 and 

(10.14) N > CG{\neN)^ 2 , 

for some sufficiently large constant C. Then (10.11) holds and moreover (3(G, N) < 
36ct/C, and so (3(G,N) is small provided that C is large enough. Choose C such 
that fi{G,N) < 1. Then c x < 2(c T + c B )/C. Choose C, so that also a < 1. 
Moreover, the bound in Theorem 6.2 is less than vkq(G + 2crG 2 (ln eiV) 1 / 2 ) < 
G2G 2 (lne7V) 1 / 2 for some constant G2, which does not depend on G, since G >> 1. 
We obtain then that the first inequality in (10.13) holds, and to satisfy the second 
inequality in (10.13), it suffices to take 

(10.15) N > eG 2 (\neN) 1 / 2 , 

for any fixed e (provided that G is large). Note that if (10.15) holds then (10.14) 
also holds. It follows that conclusion of Proposition 10.3 is true as long as inequality 
(10.15) holds for N. 
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11. Another approach 

We present now in a concise manner an approach producing an ODE similar to 
(4.9) for the determining modes Pnu, which however can also be obtained for other 
determining parameters such as nodal values, finite elements, local spatial averages 
and interpolants on fine enough spatial grids (see, e.g., [3, 21, 20, 24], for works 
related to various determining parameters). We first outline this approach for the 
case of determining modes so it can be compared to the first approach. Then we will 
state some results for other interpolants. This unified approach is fully developed 
in [12]. 

Let it(-) be a trajectory in A, and suppose w satisfies the equation 

(11.1) — — h I'Aw + B(w, w) = / — vk%uPn(w — u) , 
as 

for s > sq , with w(sq) £ H. Here [i > is a damping dimensionless parameter to 
be determined later in terms of the Grashof number G. 
The difference 6 = w — u then satisfies 

dS 
ds 

and hence 

1 d 



+ vAS + B(8, u) + B(u, 6) + B(5, S) = -vk.q(iP n 5 



2„U|2 , ,,^2„\n..X\2 



(N + iy 



< -„^\6\ 2 + -J±—\\Q N 5\\ 



Now take [i large enough so that 

(11.2) a = l-2c^G 2 + 2^>0 

and then N, large enough, so that 

< 1L3 > wrw £ r 

We then have that 



so that 



and 



~i*i a + v\\s\\ 2 < cl\5\¥\\\h ^im 2 + j^r\\ d \\ 2 

< J|5|| 2 + ^H 2 h|| 2 -^MH 2 + ^\H 2 



±\S\ 3 + au$\5\ 2 <0, 



\S(s)\ < \S(s )\ 



If the solution w(s) to equation (11.1) is defined for all ,s G [so,oo), then S(s) —> 
as s — > oo. Moreover, if w is defined, and bounded on (— oo, oo), then S(s) = for 
all set. Thus, if we first choose fi to satisfy (11.2), and then N to satisfy (11.3), 
the only bounded solution to (11.1) on (—00,00) is w = u. 
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The analog W, of the map W defined in Section 4 is now given as follows: let 
v E Cfc(R, PnV) an d let w = W(v) be the unique solution of the equation 

(11.4) ^ + vAw + B(w) = f - ukI/j,P n (w - v) 

that is bounded on R. This solution can be constructed using a Galerkin approxi- 
mation. The uniqueness of w can be proven by the same method as in the proof of 
Lemma 4.1; To obtain an a priori estimate, first note that 

^^\\w\\ 2 + v\Aw\ 2 = {f,Aw)-vKln\P N Aw\ 2 +vkIh{v,Aw) 

so 

±\\wf + 2v\Aw\* < 2(|/| + v Ko »\\v\\ x )\Aw\ , 



and we have by Lemma 3.1 



I fl 

HI <—+li\\v\\x ■ 

VK 



It is easy to show that if v(-) = Pjyu(-) for a trajectory u(-) in A, then W(v) = u. 
This suggests that we study the following analog of (4.9) 

(11.5) ^ = P N f - vAv(t) - P N B(W(v(t))) . 

The fact that for large \i the map W : v i— > W'(d) is Lipschitz is again obtained by 
a method similar to the one used in the proof of Corollary 5.2, and that the right- 
hand side of (11.5) is locally Lipschitz as a map from the space Cb(R, PnV) to itself 
is proved similarly. With these considerations the differential equation (11.5) is an 
ODE in the space Ch(R,PjvV). Note that the ODE (11.5) depends not only on N, 
but also on the parameter fi through the function W. 

As in Section 7, translation invariant solutions of (11.5) are exactly those that 
come from the NSE. For example, if v is translation invariant, then 

v(t)(s) =v(t + s)(0) =v (t + s) and w(t)(s) = w(t + s)(0) = w (t + s) . 

Moreover, since 

dv dv 

— = — at s + 1 , 
as at 

we rewrite (11.5) as 

av 

(11.6) -r + vAv + P N B{w) = P N f , 

as 

and subtract P N applied to (11.1) 
dPNW 



ds 



to obtain 



+ vAPnw + PnB(w) = PNf — vk (iPn(w — v) 



(11.7) ^+^ = 1^4(4, 

ds 

where £ = v — Pnw = Pn(v — w). Since £(s) is defined and bounded for all s G R, 
and (11.7) is linear, we have £(s) = for all s G R, provided that Kq/j, is not an 
eigenvalue of A, which can always be achieved. As a consequence (11.1) reduces to 
the NSE. 

Note that we do not have to prepare the nonlinear term as was done in (4.9) 
with the cut-off function tp in order to establish these properties. Preparation of 
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the nonlinear term is, however, needed if we want to prove the existence of absorbing 
ball for (11.5). 

One advantage that (11.5) has over (4.9) comes in the study of stationary solu- 
tions. The analog for (11.5) of the DAE in Section 10 is 

vAv + P N B{w) = P N f 

— — h uAw + B(w) = f — vk%,u,Pn(w — v) . 
as 

Note that the algebraic relation can be solved explicitly 

v = iv- 1 A- 1 P N (f-B(w)) . 

Thus in contrast with (5.3), stationary solutions of (11.5) are directly associated 
with the globally bounded solutions of an explicit Navier-Stokes type equation in 
H, namely 

(11.8) ^ + v[A + n 2 Q [iP N ]w + + kI^P n A^P n ]B{w) = [I + kIuPnA^Pn}/ . 

The study of the equations (11.5) and (11.8) for various types of determining pa- 
rameters will be presented in [12]. 

In [12] we also study a larger class of determining parameters. In particular as 
those described below. Let J be a bounded linear operator on H (approximating 
the identity) satisfying the following conditions: 

(i) \<j> — J{4>)\ < cih 2 \A(j>\, for every <f> £ D(A), where h is the size of the grid; 
moreover, h 2 is required to be at most 0((K§Glog(l+G)) _1 ), i.e. N = (kqK)~ 2 
is of the order of 0{G log(l + G)). 

(ii) JHcD(A 2 ) 

(iii) there exists a linear bounded operator K on H such that J A = K on D(A). 

Such J(u) is a determining system of parameters for (2.5). This is proven as before 
by using the equation 

(11.9) ^ + vAw + B(w, w) = f - livk 2 JJ{w) - J(u)), 
as 

where /x is in the appropriate interval. 

The map W is defined using this equation instead of (11.1). Now the equation 
which replaces the equation (11.5) is the following 

(11.10) j t v{t) + vJAW{v{t))-vKQ 2 P JH A 2 (JW{v{t))~v{t)) + JB{W{v{t))) = J} 

which is an ordinary differential equation in the ball B(0, R(rj, //)) of the space 
X = C&(M, J(H)). In (11.10) Pjh denotes the orthogonal projection of H onto the 
closure of the range of J. 

We conclude this summary by observing that given a finite rank linear operator 
J on H satisfying the condition (i), then for every sufficiently small e > there 
exists a J' satisfying all conditions (i)-(iii) and such that the norm of J — J' as an 
operator on H is less then e. 
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